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Arepeating fast radio burst associated with a 
persistent radio source 
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Thedispersive sweep of fast radio bursts (FRBs) has been used to probe the ionized 
baryon content ofthe intergalactic medium!, which is assumed to dominate the total 
extragalactic dispersion. Although the host-galaxy contributions to the dispersion 
measure appear to be small for most FRBs?, in at least one case there is evidence for an 
extreme magneto-ionic local environment** and a compact persistent radio source’. 
Here we report the detection and localization of the repeating FRB 20190520B, which 
is co-located with a compact, persistent radio source and associated with a dwarf host 
galaxy of high specific-star-formation rate at a redshift of 0.241 * 0.001. The estimated 
host-galaxy dispersion measure of approximately 90372 parsecs per cubic 
centimetre, which is nearly an order of magnitude higher than the average of FRB host 


® Check for updates 


galaxies^6, far exceeds the dispersion-measure contribution of the intergalactic 
medium. Caution is thus warranted in inferring redshifts for FRBs without accurate 
host-galaxy identifications. 


Fast radio burst (FRB) 20190520B was detected with the Five- 
hundred-meter Aperture Spherical radio Telescope (FAST)’ in drift-scan 
mode as part of the Commensal Radio Astronomy FAST Survey 
(CRAFTS)® at 1.05-1.45 GHz in 2019. Four bursts were detected during 
the initial 24-s scan. Monthly follow-up tracking observations between 
April 2020 and September 2020 detected 75 bursts in 18.5 h with a mean 
pulse dispersion measure (DM) of 1,204.7 + 4.0 pc cm”. Assuming a 
Weibull distribution of the burst waiting time, we model the FRB burst 
ratetobeR=4.5°12 h! fora fluence lower limit of 9 mJy msandaburst 
width of 1 ms, which indicates that the FRB can episodically have a high 
burst rate. Similar to other repeating FRBs, this FRB shows a complex 
frequency-time intensity structure’ with multicomponent profiles, 
sub-burst drifting and scattering (Fig. 1, Methods). No linear polariza- 
tion was detected for FRB 20190520B from the FAST observations 
(see Methods), and the rotation measure (RM) for this source has been 
detected from the higher-frequency band". The properties of 
FRB 20190520B are shown in Table 1. 

We localized FRB 20190520B with the Karl G. Jansky Very Large Array 
(VLA) using the ‘realfast’ fast transient detection system". Throughout 
the second half of 2020, we observed the source for 16 hand detected 
3,5and 1 bursts in bands centred at 1.5 GHz, 3 GHz and 5.5 GHz, respec- 
tively. We measured a burst-source position in the International 


Celestial Reference Frame (ICRF) of (right ascension (RA), declina- 
tion (dec.))(J2000) = (16 h 02 min 04.272 s, -11° 17’ 17.32") with a posi- 
tional uncertainty (1o) of (0.10", 0.08") dominated by systematic effects 
(Methods). Deep images using data from the same observing cam- 
paign revealeda persistent radio continuum counterpart at (RA, dec.) 
(J2000) = (16 h 02 min 04.261 s, -11? 17’ 17.35") with 1o position uncer- 
tainty of (0.10", 0.05") that is compact (less than 0.36"; Methods) and 
has a flux density of 202 + 8 wy averaged over a span of about 2 months 
from 30 August to 16 November 2020 at 3.0 GHz. Using the average 
flux density of each sub-band over the VLA campaign, we found that 
the radio continuum counterpart spectrum can be fit with a power-law 
spectral index of -0.41 + 0.04 (Methods). 

To identify the host galaxy, an optical image (R' band) was obtained 
using the Canada-France-Hawaii Telescope/MegaCam. Figure 2 shows 
the FRB location compared with deep optical, near-infrared (NIR) and 
radio images of the field. The optical image (R’ band) reveals the gal- 
axy J160204.31-111718.5 at the location of the FRB, the light profile 
of which peaks at about 1” southeast. Given the measured offset of 
1.3” between the FRB and the peak of the galaxy light profile, and the 
sky surface density of galaxies with this magnitude, we estimate a 
chance coincidence probability of 0.8% (Methods), indicating that 
J160204.31-111718.5 is the host galaxy of FRB 20190520B. The NIR 
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Fig.1|BurstsfromFRB 20190520B, shown as dynamic spectra and burst 
profiles. Top: frequency-integrated burst normalized intensities that were 
detected during FAST observations. Thesesixbursts are chosen from different 
observation epochs. The burst labels ‘PXX’ are in order of arrival time and the 
corresponding burst properties are shownin Supplementary Table 1. Bottom: 
de-dispersed dynamic spectra of the bursts, clearly showing their band-limited 


image (J band, 1.153-1.354 um) obtained with Subaru and Multi-Object 
Infrared Camera and Spectrograph (MOIRCS)” most likely shows the 
stellar continuum emission of the host galaxy with an AB magnitude of 
22.1+ O.1inthe] band, with FRB 20190520B and the radio continuum 
counterpart on the galaxy periphery. 

We obtained an optical spectrum at the location of the FRB with 
the Double Spectrograph on the Palomar 200-inch Hale Telescope, 
which reveals the redshift of the putative host to be z= 0.241 + 0.001 
based ona detection of strong Ho, [O III] 4,859 À and [O III] 5,007 À lines 
(Methods). Afollow-up observation with the Low Resolution Imaging 
Spectrometer (LRIS) at the Keck I Telescope covering both the FRB 
location and the nearby Subaru J-band source along the extended 
R’-band structure indicates that the R’-band structure is dominated 
by the [O III] emission at the same redshift of z = 0.241. The Ha lumi- 
nosity Li, = 7.4 € 0.2 x 10*° erg s ! after extinction correction implies 
astar-formation rate of about 0.41 M, yr™. On the basis of the J-band 
magnitude, we estimate the stellar mass of the host galaxy to be about 
6 x 10° M,. Thus, we characterize J160204.31-111718.5 as a dwarf galaxy 
with a relatively high star-formation rate for its stellar mass compared 
with local Sloan Digital Sky Survey (SDSS) galaxies". At the luminosity 
distance of 1,218 Mpc implied by the redshift, the radio continuum 
counterpart has a spectral luminosity of L44,, = 3 x 10? erg s” Hz". 

The redshift of an FRB source is often estimated from the 
DM attributed to the intergalactic medium (IGM), DMigy = 
DMerp - DMpost - DMyw Where MW is Milky Way. Theoretical calcula- 
tions? and observations! have independently estimated the IGM 
contribution to FRB DM as a function of redshift (Fig. 3). Fora DM 
contribution from the Milky Way of 113 + 17 pc cm? (including a +40% 
range for the NE2001 model estimate” of the disk contribution and 
auniform halo distribution from 25 pc cm?to 80 pc cm?) combined 
with an assumed host-galaxy DM of only 50 pc cm”, the implied red- 
shift range for FRB 20190520B is z = 2.2 toz = 0.9 for baryon fractions 

fiom Of 0.6 to 1 for the ionized IGM, much larger than the measured 
value (Methods). 

Instead, using DMigy(z= 0.241) «195*I?pccm ? (68% interval to 
account for the cosmic variance of DM,,y and 0.85 for the ionized 
baryon fraction") combined with the Milky Way contribution, a very 
large host DM (disk + circumgalaxy + local to source contribution) of 
DMpost = 90372 pccm ? is inferred from the posterior distribution. 
Over a broader range of ionized fractions from 0.6 to 1, DM, extends 
from 1,020 pc cm? to 745 pc cm”. 
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nature. The DMsare taken from Supplementary Table1. The colour map is 
linearly scaled and darker patches represent higher intensities. The bad 
frequency channels are masked and labelled using red patches onthe left. The 
time and frequency resolutions for the plots are downsampled to 0.786 ms and 
3.91 MHz, respectively. 


In addition to the low chance coincidence probability, the measured 
DM and scattering properties of FRB 20190520B render it unlikely 
that J160204.31-111718.5 is a foreground galaxy with the true FRB 
host much farther in the background. The observed Ha emission 
implies a DM contribution from warm gas in the galaxy ranging from 
230 pc cm? to 650 pc cm” (observer frame) for temperatures from 
0.5x10*Kto 5x 10* K (Methods). The scattering contribution of a 
foreground galaxy depends not only on its DM contribution but also 
onthe geometric leverage effect, which will increase the scattering by 
several orders of magnitude relative to scattering in the host galaxy. 
If J160204.31-111718.5 were a foreground galaxy, we estimate a fidu- 
cial range of scattering times from 0.6 s to 20 s at 1.25 GHz, orders of 
magnitude larger than the observed mean scattering time of 10 + 2 ms. 
Although this estimate depends on parameters such as the ionized 
gas temperature, path length through the galaxy and the degree of 
turbulence in the gas, which are not well constrained, the observed 
scattering time is probably too small for the proposed host galaxy to 
lieintheforeground (Methods). 

The large DM,,, inferred for FRB 20190520B demonstrates that the 
distribution of DM,,,. values for the FRB population has a long tail, which 
may add considerable varianceto estimates for the IGM contribution. 
It is conceivable that the DM,,, distribution may differ for repeating 
and non-repeating FRBs, which could make non-repeating FRB DMs 
more accurate proxies for redshift. To understand this further, a larger 
sample of precisely localized FRBs with measured host-galaxy redshifts 
is needed to statistically characterize host galaxies, their circumgalactic 
media and near-source environments along with the IGM”. Accordingly, 
searches need to accommodate large values of DM, as part of the DM 
budget for FRB sources. Large DM, contributions may also imply large 
scattering, which can also reduce search sensitivity. 

The co-located radio continuum counterpart, the star-forming 
dwarf host galaxy and the high repetition rate make FRB 20190520B 
a clear analogue to FRB 20121102A, the first known repeating FRB” 
and the first to be identified with a compact, luminous persistent 
radio source (PRS)?? (I, (4,7 3 x 10? erg s ! Hz”). Another repeating 
source, FRB 20201124A, was also associated with a radio continuum 
counterpart” ?*: however, optical spectroscopy and radio interfero- 
metric measurements demonstrated that the persistent radio emission 
was spatially extended and consistent with an origin in star forma- 
tion in the host galaxy”. By contrast, the continuum counterpart to 
FRB 20190520B appears to not be from star formation because its 


Table 1| Properties of FRB 20190520B 


Burst parameters 


Right ascension (J2000) 


16 h O02 min 04.272 s (0.1") 


Declination (2000) 


-11° 17 17:32" (0.08") 


Galactic coordinates (l, b) 


359.67°, 29.91? 


Number of detections? 88 

Mean total DM (pc cm ?) 12047 + 4.0 
Mean burst width (ms) 13.5 € 1.25 
Mean scattering timescale(ms)at1.25GHz 10+2° 
Measured fluence (Jy ms) 0.03 to 0.33 
DMuw,disk- DMyw, halo (pc cm 2 60, 50 
DMhost (pc cm?) 9032 
Luminosity distance (Mpc) 1,218 
Isotropic equivalent energy (10?! erg) 3.6 to 40 


Persistent radio source 


Right ascension (J2000) 


16 h O2 min 04.261 s (0.1") 


Declination (2000) 


-11° 171735" (0.055 


Flux density at 3.0 GHz (uJy) 


2028? 


Luminosity at 3.0 GHz (erg s” Hz") 


3x10? 


Size (3 GHz) 


«0.36" («1.4 kpc) 


Host galaxy 

Redshift (z) 0.241 + 0.001 
Stellar mass (M,)° ~6 x 108 

Ha luminosity (erg s”) 7.4 x 10%? + 0.2 x 10? 
Star-formation rate (M, yr’) ~0.41 


*Including the FAST and VLA observations. 


PThese are the mean values from FAST bursts fitted with a Gaussian pulse convolved with a 


one-sided exponential (Methods). 


“The MW electron density model from NE2001 and YMW16. 
dAveraged over a span of about 2 months from 30 August to 16 November 2020. 


*Based on the J-band magnitude. 
‘Based on the Ha luminosity. 


luminosity would imply a star-formation rate of about 10 M, yr !, a 
factor of 25 larger than that measured for the host galaxy anda factor 
of five larger than the highest observed star-formation rate for galaxies 
of this mass". Given its extreme luminosity, unresolved structure in 


a b 


VLA observations and offset from the peak optical emissionin the host 
galaxy, we conclude that the radio continuum counterpart is a compact 
source (less than 1.4 kpc) physically connected to FRB 20190520B, 
a PRS like that associated with FRB 20121102A (which has a very 
long baseline interferometry (VLBI) confirmed size of less than 0.7 pc 
(ref. ?9)). 

More than a dozen FRBs were localized before FRB 20190520B, 
including five repeating sources”, but only FRB 20121102A had been 
associated with a compact PRS. FRB 20121102A also demonstrates a 
sporadically large burst rate” (for example, a peak burst rate of 122 h ?), 
a substantial RM that varies over both short (burst to burst) and long 
(year) timescales** and DM,os as large as about 300 pc cm”, suggest- 
ing that burst activity may be correlated with both relativistic plasma 
emitting synchrotron radiation and the presence of thermal plasma in 
thelocal FRB environment. In addition to a PRS, FRB 20190520B shows 
a DM, that is almost three times larger than that of FRB 20121102A, 
and may have a comparably large RM (Methods). If a sizeable frac- 
tion of DM, is from thermal plasma in the circum-source medium, 
then perhaps the presence of a PRS and a dynamic magneto-ionic 
environment are correlated with FRB formation, and repeating FRBs 
like FRB 20121102A and FRB 20190520B are younger sources that still 
reside in their complex natal environments. 

However, until further VLBI observation constraints become avail- 
able, it is unclear how much of the large DM,,, for FRB 20190520B is 
attributable to ionized gas in the circum-source medium versus the 
host galaxy's interstellar medium. Moreover, other repeating FRBs 
have very deep limits on PRS counterparts????, which complicates the 
connection between burst and magneto-ionic activity and PRS lumi- 
nosity. For example, FRB 20200120E has an upper limit on persistent 
radio luminosity of, sc, < 3.1x 10? erg s! Hz ! (ref. ?).Itis associated 
with an old stellar population and has a modest DMpost < 50pccm ^, 
which is in sharp contrast to FRB 20121102A and FRB 20190520B. One 
possibility is that there are multiple kinds of sources that can emit FRBs, 
a point that has been argued based on burst rate and phenomenol- 
ogy???! Alternatively, the source properties may evolve as the source 
ages, the PRS fades, the event rate drops and the surrounding plasma 
dissipates”. The similarity of FRB 20190520B to FRB 20121102A sug- 
gests a potential connection between burst activity, the presence of a 
PRS and DM,,,. for at least some FRBs. 

Various methods have been used to argue that repeating and 
non-repeating FRBs comprise different subclasses????, either at dif- 
ferent evolutionary phases or entirely different physical scenarios. 
Although the observed burst repetition and morphology can be time 


c 


Fig.2|Optical, infrared and radio images ofthe field of FRB 20190520B. 

In each case, the box is 40"in size and the 2” crosshairs indicate the best FRB 
position at RA - 16 h 02 min 04.272 s, dec. = -11° 17’ 17.32" (J2000). North is up 
andeastistotheleft. a, An optical R-bandimage obtained by CFHT MegaCam 
covers 5,427-7,041 À, including redshifted Hf 4,861 A, [O III]4.959 À and 

[O 111]5,007 A emission lines from the host galaxy. The bright streak in the lower 
left is an artefact caused by a bright star outside the field of view. The inset 


shows the host galaxy ina5” region centred onthe FRB position, as indicated by 
the yellow circle (0.1" radius, corresponding to thelo position uncertainty). 

b, Theinfrared]-band image by Subaru/MOIRCS shows emission only at the 
location of the peak of the optical light profile of the host galaxy. The inset is a 
5” region matching the inset ina.c, The radio VLA image (2-4 GHz) shows a 
compact persistent source at the FRB location. The synthesized beamis shown 
as anellipse of size (0.92” x 0.47”) inthe left corner. 
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Fig. 3 | Galactic and extragalactic contributions tothe DM observed for 15 
FRBs. Theses 15 FRBs have firm host-galaxy associations and redshifts. 
Sources localized on initial detection (‘one off’ bursts) and repeating FRBs are 
identified separately. Galactic disk contributions are estimated from NE2001/* 
with +20% uncertainties, along with an additional halo contribution of 

25-80 pc cm” (full range). The red error bars on each extragalactic DM 
estimate represent aconservative full range uncertainty derived by adding the 


dependent owing to various mechanisms”, PRS emission and DM, ost 
may reflect more persistent aspects of the FRB environment and 
thus may be more reliable tracers of any putative subclasses. Further 
progress will result from better characterization of the full range of 
host-galaxy DMs, which will also help mitigate biases in the DM-red- 
shift relation’. 
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Methods 


Observations 

FAST. CRAFTS is a multi-purpose drift scan survey conducted with 
FAST using a 19-beam receiver operating at 1.05-1.45 GHz, deployed 
in May 2018 and conducting blind FRB searches using multiple pipe- 
lines?*. FRB 20190520B was discovered on 16 November 2019 in archived 
CRAFTS data, which are in 1-bit filterbank format with 196-pis time reso- 
lution and 0.122-MHz frequency resolution. In this first discovery obser- 
vation, 3 bursts were detected in 10 s, and another burst was detected 
20 slater. These 4 bursts from the drift scan survey gave a preliminary 
location for the source within an approximately 5-arcmin-diameter re- 
gion. Taking the pointing location from FAST, 2 follow-up observations 
were performed with FAST on 25 April and 22 May in 2020 with 19-beam 
mode in which 15 bursts were detected. A monthly observation cam- 
paign was then conducted by FAST using the approximately 100-mas 
localization from the VLA (next section). After some regular telescope 
maintenance, 10 observations were performed spanning from 30 July 
2020 to 19 September 2020 in which 60 more bursts were detected. 

Bursts were detected from FRB 20190520B in each FAST monitoring 
observation. We list the properties of those bursts detected by FAST in 
Supplementary Table 1. The burst arrival time is in modified Julian date 
(MJD) format and has been transformed to the arrival time at the Solar 
System barycentre at 1.5 GHz with the DM values from Supplementary 
Table 1. The observed DM,;, is measured using the method from ref. ? 
and the code from ref. ?. We use the DM value that maximizes the struc- 
ture of the highest signal-to-noise ratio (S/N) burst on a given date, 
and report this same DM for all other bursts detected on that date. 
Apparent epoch-to-epoch DM variations may simply be the result of 
a combination of variable effects, including intrinsic pulse structure, 
time-frequency drift, flux distribution in frequency and scattering 
(which affects bursts chromatically). 

The FAST data were searched using a Heimdall-based pipeline****. For 
the FRB blind search in the 19-beam drift scan survey, the polarizations 
were added and only the total intensity (Stokes I) was recorded. The 
trial DM range was from 100 pc cm ?to 5,000 pc cm? and we matched 
the pulse width using a boxcar search up to 200 ms. The candidates 
with S/N > 7 and present in fewer than 4 adjacent beams were manually 
examined for further inspection. After we identified FRB 20190520B 
as anew FRB source, the follow-up burst search was done with a nar- 
rower DM range (100-2,000 pc cm). Following the VLA localization, 
the tracking observations only recorded data from the central beam 
but with all the Stokes parameters and higher time resolution (about 
50 us). Candidates with S/N > 7 correspond to a fluence threshold of 
9 mJy ms fora burst width of lms. 

The pulse widths were estimated by a Gaussian profile fit if the burst 
showed no evidence of scattering (See section on scattering below). 
A subpulse is recognized if the bridge between two closely spaced 
peaks drops more than 5o below the higher peak. If the burst shows a 
scattering tail, the pulse width is derived from a combined fit for the 
Gaussian width and scattering time. We roughly estimate the bandwidth 
of each burst by dividing the whole bandpass into 50-MHz subbands 
and identifying the subbands containing burst emission. The burst 
fluences were determined using the bandwidth, temporal width and 
amplitude of each burst. 


VLA. Following the FAST detection of bursts from FRB 20190520B, VLA 
observations were performed from July-November 2020 with direc- 
tor's discretionary time (DDT) project 20A-557 at the most reliable 
position determined using FAST (good to 5 arcmin). Most of the ob- 
servations were in the B array configuration (with a maximum base- 
line of 11.1 km), with the exceptions that the array configuration was 
BnA on MJD 59,161 and BnA->A on MJD 59,167 and MJD 59,169. In total, 
11.4 hwere spent onthe source at multiple bands. The total bandwidths 
for the L (1.5 GHz), S (3 GHz) and C (5.5 GHz) bands were 1,024 MHz, 


2,022 MHz and 2,022 MHz, respectively, with 1,024 channels, corre- 
sponding to channel bandwidths of approximately 1 MHz, 2 MHz and 
2 MHz, respectively. 

The details of the observations are given in Extended Data Table 1. 
The telescopes were pointed at the field centred at (RA, dec.) 
(J2000) = (16 h 02 min O1s, -11? 17’ 28"). We used the ‘realfast’ search 
system at the VLA to search for bursts from FRB 20190520B in our VLA 
observations. Observations on MJD 59,169 contributed to the imaging 
at S band, but the ‘realfast’ system was not run on this day owing toa 
system error. 

The ‘realfast’ search system has been described in detail in refs. 
but here we discuss it briefly. Using a commensal correlator mode, 
the visibilites are sampled with 10-ms resolution and distributed 
to the ‘realfast’ graphics-processing-unit cluster to search for tran- 
sients. The search pipeline rfpipe** then applies the online calibration, 
de-disperses the visibilities and forms images at a set of different 
temporal widths. The 80 fluence limit of a 10-ms image is 0.29 Jy ms, 
0.18 Jy ms and 0.13 Jy ms for the L, S and C bands, respectively. Can- 
didates with an image S/N greater than the threshold triggered the 
recording of fast sampled visibilities in windows ranging from 2 s to 
5scentred on each candidate. For each candidate, time-frequency 
data corresponding to the maximum pixel in the image are then 
extracted and post-processed. The candidate is then classified by 
fetch, a graphics-processing-unit-based convolutional neural net- 
work”. 

Promising candidates selected visually by the ‘realfast’ team go 
through an offline analysis to refine the candidate parameters and to 
improve its detection significance, if possible. Several methods are 
tried within this analysis: offline search for the transient using a finer 
DM grid, varying radio frequency interference (RFI) flagging thresholds, 
changing image size (as the real-time search is done on non-optimal 
image sizes owing to computational limitations), subband search and 
so on (see section 2.4 of ref." for more details). The pixel sizes of images 
formed by the real-time system for the L, S and C bands were 0.9”, 0.48" 
and 0.27”, respectively, with an image size of 2,048 x 2,048 pixels cor- 
responding to field sizes of 0.5°, 0.27° and 0.15°, respectively. During 
refinement, we searched the data with smaller pixel sizes: 0.5” x 0.75” 
for the L band, 0.38” x 0.28” for the S band and 0.27” x 0.18” for the C 
band. The significance of an astrophysical transient improves when 
the data is de-dispersed ata DM closer to the true DM of the candidate 
(see section 2 of ref. *°). Therefore, we re-ran the search witha finer DM 
grid at 0.1% fractional sensitivity loss, compared with 5% used by the 
real-time system. Noise-like events or RFI are sensitive to RFI flagging 
and image gridding parameters, and so they cannot be reproduced 
on refinement and are discarded. Sometimes, the transient signal is 
present only ina fraction of the whole frequency band, so we also re-run 
our search only using the relevant frequencies, which further improves 
the detection significance. We applied these techniques on all the can- 
didates selected from these observations to refine their parameters. 
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Localization of bursts 

Calibration. The ‘realfast’ system makes several assumptions during 
calibration and imaging to improve the computational efficiency for 
the real-time imaging and search. Moreover, the point spread function 
(PSF) of the interferometer is not deconvolved from the image, so the 
real-time system forms ‘dirty’ images. The PSF shape makes the images 
more difficult to visually interpret and model. 

To address these issues, we used a top-hat function to select the 
raw, de-dispersed, fast-sampled visibilities containing only the burst 
signal to re-image the burst data using the Common Astronomy 
Software Application package (CASA v5.6.2-3)". Observations of 
3C 286 (acquired before the FRB observations) were used to calibrate 
the flux density scale, bandpass and delays. Complex gain fluctua- 
tions over time were calibrated with observations of the calibrator 
J1558-1409. We performed phase-reference switching at intervals 


of 16 min, 13 min and 12 min for the L, S and C bands, respectively, 
consistent with the nominal phase-coherence timescales for the VLA, 
so that after calibration any short-timescale phase variations are 
negligible. Therefore, the systematic errors for the burst positions 
on short timescales are ofthe same magnitude as those for the deep 
imaging, as discussed below. 


Determining properties of individual bursts. After calibration, the 
CASA task tclean was used to generate an image for each burst and 
estimate the S/N. Most of the bursts are spectrally confined, so for each 
burst we select the spectral window range that produces the highest 
image S/N. We then use the CASA task imfit around the FRB position to 
fit an elliptical Gaussian to the source in the radio image and measure 
the centroid location, peak flux density and loimage-plane uncertain- 
ties. Hereafter, we refer to these image-plane uncertainties as statistical 
errors on burst positions. 

The average burst positions are calculated for each frequency band 
separately, by weighting each detection by the inverse of the position 
fit errors (statistical errors) reported by CASA. Statistical errors are 
inversely proportional to the S/N, and therefore high significance detec- 
tion is expected to have smaller fit errors. The total positional error at 
each frequency band is obtained by adding this statistical error in 
quadrature with the systematic error in that band (determined using 
the deep radio images described below). The total positional error (RA 
error, dec. error) is (0.25”, 0.32”) for the L band, (0.28”, 0.17”) for the S 
band and (0.12”, 0.09”) for the C band. In the final analysis, the errors 
are dominated by systematic uncertainties at each frequency band 
(Extended Data Tables 2, 3). The final burst position was obtained by 
taking a weighted average of the burst positions at each of the three 
frequency bands, using the inverse of the total error at each frequency 
band as the weight. The best estimate of the burst position is 
RA -16 h 02 min 04.272 s, dec. = -11? 17’ 17.32" (J2000). We estimate 
the error on this position to be 0.10” and 0.08”. We calculated the 
reduced chi-square value of individual burst positions with respect to 
this best estimate of the FRB position. The reduced chi-square was 
given by y? = (1/8) Y; (0;- Oves) /02, Where 6, refers to the coordinates 
ofa burst, 0,.., refers to the best estimate of the FRB position, o; is the 
total error ona burst (statistical and systematic errors added in quad- 
rature) and the sum is over all the 9 localized bursts. We obtained a 
reduced chi-square value of 0.53 and 0.77 for RA and dec. respectively. 

We confirmed that there were no short-timescale phase errors by 
performing intermediate-timescale imaging of the continuum data. 
We imaged segments of 5-30 s, such that at least 1-3 sources could be 
detected at each frequency, and inspected the position variations of 
those sources over time. At all frequencies, the positions were stable 
and consistent with the statistical uncertainty (radiometer noise). 
Offsets of the sources were within the range of the systematics quoted 
for the PRS in the next section. 

Owing to the lowtime resolution of the VLA data (10 ms), we do not 
perform any sophisticated modelling of the burst properties. The prop- 
erties of the VLA bursts, including the flux densities, were obtained 
using CASA and the ‘realfast’ system, and are given in Extended Data 
Table 4. Here the reported widths should be considered as upper limits. 
The frequency extent of the burst signal was visually determined and 
reported in the last column. We also could not estimate a structure 
maximizing DM for the VLA bursts because of the coarse time 
resolution. The DMs reported in this table are the values that 
maximize the S/N. This also explains the apparent variability seen in 
the DM values. 


Persistent radio source 

Deep radio images and PRS. The VLA campaign obtained two epochs 
at 1.5 GHz and six epochs at 3 GHz and 5.5 GHz, which resulted in an 
on-source time of about 3 h for 1.5 GHz, and about 4 h for both 3 GHz 
and 5.5 GHz. 


Along with the ‘realfast’ output, the VLA visibilities with 3s or 5s 
sampling times were saved and analysed to search for persistent radio 
emission. This is done in parallel with saving the high-time-resolution 
(10 ms) data around the burst that was used in the previous section. We 
used the same data reduction, flagging and calibration approach that 
was used for burst localization. We then performed further flagging 
onthe target and then subsequently imaged its Stokes I data using the 
CASA deconvolution algorithm tclean. To balance sensitivity while 
reducing sidelobes from a nearby bright source, we imaged witha 
Briggs weighing scheme (robust=0). In addition, self-calibrations were 
performed for all observations to correct considerable artefacts from 
the close-by bright sources in the field. We made use of the CASA task 
imfitto measure source flux densities by fitting an elliptical Gaussian 
model in the image plane. 

We stacked observations at each central frequency in the u-v plane 
and then imaged the Stokes I intensity, resulting in deep images at 
1.5 GHz, 3 GHz and 5.5 GHz with root-mean-square (r.m.s.) noise of 
9.0 Wy per beam, 4.5 Wy per beam and 3.0 py per beam, respectively. 
At 1.5 GHz,3 GHz and 5.5 GHz, the upper limits of the deconvolved sizes 
of the PRS are as large as 1.4” x 0.89", 0.51" x 0.14” and 0.36” x 0.1”, 
respectively. Thus, a conservative upper limit of the size of the PRS 
from VLA observations is 0.36", which corresponds to 1.4 kpc at the 
angular diameter distance of 809 Mpc. The obtained positions for the 
PRS are shown in Extended Data Table 2 and the positions of the bursts 
and PRS are shown in Extended Data Table 1. The systematic offsets on 
these positions are estimated in the next section. 


Systematic offsets. To determine the systematic errors on the coordi- 
nates of the PRS that we determined from the deep images, we ran the 
PyBDSF (https://www.astron.nl/citt/pybdsf/index.html) package to 
extract radio sources from the deep images. We then cross-matched the 
detected point sources in the deep images with the sources listed in the 
optical PanSTARRS survey DRI?. We identified the radio point sources 
using the following criteria: (1) the peak intensity (Jy per beam) ofa 
source should be 0.7, 0.5 and 0.5 times higher than its integrated flux 
(Jy) forthe1.5-GHz, 3-GHz and 5.5-GHz images, respectively; (2) the S/N 
(peakintensity/local r.m.s. noise) ofa source should be greater than 5. 

In total, we detected 375, 113 and 43 sources in the 1.5-GHz, 3-GHz 
and 5.5-GHz deep images, respectively. The selected sources were also 
visually checked to make sure that they are ‘point like’ sources in the 
deep images. We searched for matching optical sources within radii of 
0.5”, 1.0" and 2.0". Going from 1.0" to 2.0”, the additional cross-matched 
sources at each band are consistent with chance coincidence. There- 
fore, we adopt a 1.0" cross-match radius, finding 136, 31and 9 sources 
with optical counterparts in the PanSTARRS DRI catalogue at 1.5 GHz, 
3 GHz and 5.5 GHz, respectively. These cross-matched sources were 
used to determine the astrometry for FRB 20190520B and its PRS. By 
subtracting PanSTARRS coordinates from VLA coordinates, we estimate 
the systematic offsets in the 1.5-GHz, 3-GHz and 5.5-GHz positions, as 
listed in Extended Data Table 2. The systematic offsets are consistent 
with zero mean, and their uncertainties dominate the uncertainties 
ofthe PRS position. 


Determining the position of the PRS. To determine the position of 
the PRS in the three frequency bands, we followed a procedure simi- 
lar to what was used to determine the burst positions (see ‘Localiza- 
tion of bursts"). The best estimate for the PRS position is (RA, dec.) 
(J2000) = (16 h 02 min 04.261 s, -11? 17’ 17.35”). The error on this posi- 
tion is estimated to be 0.10” and 0.05” for RA and dec., respectively. 


Variability and spectrum of the PRS. The flux density of the source 
measured at each epoch is shown in Extended Data Table 1. The meas- 
ured flux densities show variations that are mostly consistent with 
measurement errors, but there are roughly 2o variations in the 2.5-GHz 
subband that, if real, could be refractive interstellar scintillation or 
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intrinsic variations, or both. To study the spectrum of the PRS, we split 
each of the observations into two 0.5 GHz/1 GHz subbands. Then we 
measured the average flux density at each of the subbands over the 
campaign. The multiband data were fit with a power-law model (S, = v^, 
where S, is the observed flux density at the frequency v, and a is the 
spectral index), yielding an average spectral index for the PRS of 
-0.41 + 0.04 (Extended Data Fig. 2). 


Chance coincidence association of the PRS. In the 1.5-GHz deep 
image, we detected 8 ‘point like’ sources, including the PRS, within 
5 arcmin of the phase centre and with a flux density no less than that 
of the PRS (at 260 uJy) based on our point-source selection criteria de- 
scribed earlier. There is an additional bright source with a flux density 
of a few tens of millijansky in the region, but it was not classified as a 
‘point like’ source based on our criteria. 

To estimate the chance coincidence probability of the PRS with the 
bursts, we compared the solid angle corresponding to the uncertainty 
of burst localization and the average solid angle occupied by each of 
the eight ‘point like’ sources in the FRB field of view. The solid angle 
corresponding to each of the 8 ‘point like’ sources is roughly estimated 
as Sou, = 11(5/60)/8 steradians (Sr). The offset between the average 
position of the nine bursts and the position of the PRS at 5.5 GHz, which 
is best constrained when taking both statistical and systematic errors 
into account, is about 0.06 arcsec. This, along with a statistical error 
of 0.01 arcsec and a conservative estimate of the systematic error of 
0.12 arcsec, can be used to estimate the offset between the PRS and the 
FRB position. We conservatively estimate this offset to be 0.19 arcsec. 
The solid angle corresponding to the offset therefore can be estimated 
as Soc 11((0.19/60/60)"/8 Sr. The ratio between 5,4... and Ssource gives 
the chance of coincident association of the PRS with the FRB position 
to be about 3 x 10°. 


Galaxy photometry and redshift determination. The deep R’-band 
(5,427-7,041 A) images obtained by CFHT/MegaCam were stacked from 
archival observation data taken in 2014-2015 by the Canada France 
Hawaii Telescope (CFHT) archival pipeline MEGAPIPE, with a total of 
about 3.6 hon the field. The level 3 (flux calibrated) images were re- 
trieved for our analysis. 

NIR J-band images of the FRB-20190520B field were taken under 
a relatively poor seeing condition (about 1.3”) through a Subaru 
target-of-opportunity observation on 5 August 2020. A total of 1.4h 
of observations were used for the final combined J-band image shown 
in Fig. 2. AJ-band source of 22.07 + 0.14 mag (AB) was detected at about 
1” southeast of the burst location, and is possibly stellar emission from 
the host galaxy. A faint source 2.5” north has 22.87 + 0.26 mag in the] 
band. Neither of these two sources were detected in the Ks-band image, 
with a 5o limit of 21.74 mag (1.1 h). 

An optical spectrum was obtained with the Double Spectrograph 
(DBSP) on the Palomar 200-inch telescope on 24 July 2020 using a 
1" slit width. This observation was executed before the CFHT/Meg- 
aCam archival data on FRB 20190520B field were found, and only 
PanSTARRS images were used for observation planning. The slit of 
DBSP was set to cover the PRS emission at RA = 16 h 02 min 04.27 s, 
dec. = -11° 17' 17.5" detected by the VLA in L-band on 22 July 2020. No 
clear optical counterpart was detected in any of the five band images 
of PanSTARRS from DRI”. The slit was guided by the nearby M-star at 
RA =16 h 02 min 04.48 s, dec. = -11? 17' 19.1" as reported in the PanS- 
TARRS DRI catalogue, with i= 20.4 mag, and 3.4" due east of the PRS 
coordinates. Theslit was setto a position angle of 108.5?, ensuring that 
both the PRS and the M-star fell within the slit. The observations with 
2x900sexposures werecarried outunder photometric sky conditions 
and subarcsecond seeing. The two-dimensional spectrum was gener- 
ated using the image reduction and analysis facility (IRAF), including 
bias removal, flat-fielding and reduction of other instrumental effects. 
The one-dimensional spectrum was extracted from a 1.5" window. 


The standard star BD+284211 was used for telluric correction and flux 
calibration. The DBSP one-dimensional spectrum is shown in Extended 
Data Fig. 3. The flux scale of the spectrum does not include the slit loss 
and registration error of PanSTARRS coordinates of the M-type star. 
The [O III] 5,007 A line and the Ha line are both well detected (>50). 
The two emission lines are narrow, each with a full-width at 
half-maximum of about 10 A. The redshift derived from these two 
spectral lines is z= 0.241 + 0.001. 

A follow-up Keck LRIS spectroscopic observation was carried out on 
25 August 2020 under reasonable weather and seeing conditions (1.1"). 
The 1.5” slit was set at a position angle of 160? to the extended optical 
emission seen in the MegaCam R'-band image around the 
FRB 20190520B location. A total exposure of 3,600 s was obtained. 
The emission lines Ha, Hf, [O III] 4,859 A, and [O III] 5,007 A are well 
detected, indicating that the extended R’-band structure has the same 
redshift of z= 0.241. After the spectrum was corrected for Galactic 
foreground extinction, the line flux ratio between Ha and HB was used 
to estimate an extinction of A, = 0.80 mag, assuming a case-B line 
ratio of 2.86, yielding an Ha flux Fuya = (42.0 € 0.5) x 10 “ergs ‘cm 2 
Theextinction-corrected Ha luminosity is Lya = 7.4 0.2 x 10*°ergs 7. 


Chance association probability of the host galaxy. We use the ap- 
proach of ref. to estimatethe chance coincidence probability of the 
association between the galaxy (J160204.31-111718.5) and 
FRB 20190520B. Assuming a uniform surface distribution of galaxies, 
the probability of chance coincidence follows a Poisson distribution; 
thatis, P- 1- exp( - nj), where n;is the mean number density of galax- 
ies brighter than a specified R’-band magnitude in a circle of given ra- 
dius, determined by the half light radius of the galaxy and the burst 
localization error region. This number density is estimated using the 
results from ref. ^*. From previous discussion, we assume the localiza- 
tion error onthe FRB position to be about 0.1”, and the size of the host 
galaxy to be 0.5”. The R’-band magnitude of the possible host is difficult 
to estimate asit is significantly affected by the emission lines. Thus, we 
conservatively estimate the R’-band magnitude to be <23.3” mag as- 
suming a flat spectral energy distribution in vL, between the R’ band 
and theJ band, where v and L, are frequency and luminosity, respec- 
tively. The probability of a chance association between FRB 20190520B 
and J160204.31-111718.5 is estimated to be about 0.8%, independent 
of constraints from the association with a rare PRS or the observed FRB 
scattering. The next nearest galaxy is 6.5" away from the FRB and has a 
chance coincidence probability of >20%. 


FAST burst sample analysis 

Repetition rate. FRB 20190520B was found to be regularly active with 
two or more bursts detected in each monitoring observation with FAST. 
Owing tothe possible clustering behaviour of FRB emission, a Weibull 
distribution was used for the waiting time ó (the separation between 
adjacent detected bursts)? 
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where kis the shape parameter and A is the scale parameter of the 
distribution. The Weibull distribution reduces to a Poisson distribu- 
tion when k=1. The mean waiting time, F(6) - AI (1- 1/k), where I is 
the gamma function, implies a burst rate r=1/F(6). Extended Data 
Fig. 4 shows the inferred parameter distributions obtained with a 
Markov chain Monte Carlo. We find the burst rate of FRB 20190520B 
is r= 4.512 h! with shape parameters k = 0.3770-04 when all 79 bursts 
above 9 mJy ms are used (left panel in Extended Data Fig. 4). Waiting 
times shorter than1 sinclude substructure in individual bursts, so we 


also analyse only longer wait times 6>1s, yielding r=5.3"t} h" 


with shape parameters k= 0.76*00? (right panel in Extended Data 


Fig. 4). 


Short- and long-timescale periodicity search. For a range of trial 
periods (P) and period derivatives (P), thearrival times ofthe 79 FAST 
bursts were folded to calculate the phase of each burst for a given trial 
period. For the short-timescale periodicity search, the arrival times 
were folded at Pbetween 1 ms and 1,000 s and P between 10 ? s s! and 
1s s”. For the long-timescale periodicity search, the trial range of pe- 
riod Pwasfrom2dto365 d and the range for period derivative (to ac- 
count for spindown) was P from 10 "? d d! to 1 d d”. The longest con- 
tiguous phase segment without bursts was calculated for each fold 
trial as a signature for possible periodicity. 

For the short-timescale periodicity search, the bursts spread to a 
phase window larger than 60% of one period, which indicates no peri- 
odicity pattern in the trial range. For the long-timescale periodicity 
search, the longest phase segment without bursts showed two maxima 
about 0.6 near 67 d and 169 d. These two maxima are reproduced by 
folding the MJDs of all observation sessions, which indicates that the 
maxima reflect the sampling window rather than any burst periodicity. 
Thus, no long or short period of FRB 20190520B was detected. 


Energy distribution. A 1-K equivalent noise calibration signal was in- 
jected before each observation session to obtain a high-quality flux 
density and energy calibration measurement for each detected burst. 
The injected noise calibration signal was used to scale the data into 
temperature units, yielding anearly constant r.m.s. radiometer noise 
to within 6%. We then converted from units of kelvin to jansky using 
the zenith-angle-dependent gain curve, provided by the observatory 
through quasar measurements“. For most observations, the zenith an- 
gleis less than 20°, which corresponds toa stable gain of about 16 K Jy”. 

Assuming bursts have spectral index of about 0, we calculated their 
isotropic equivalent burst energies, £, following equation (9) of ref. ": 
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where F,=S, x W,,isthe specific fluence in units of erg cm? Hz orJy ms, 
S, is the peak flux density, which has been calibrated with the noise level 
of the baseline, giving the flux measurement for each pulse at acentral 
frequency of v,=1.25 GHz, W,, is the equivalent burst duration, and the 
luminosity distance D, = 1,218 Mpc corresponds to a redshift z = 0.241 
for the source of FRB 20190520B, using standard cosmological param- 
eters”. The fluence-width distribution at 1.25 GHz for FRB 20190520B 
bursts can be seen in Extended Data Fig. 5. The histogram of burst ener- 
gies (Extended Data Fig. 6) exhibits abump that we fit with alog-normal 
function. 


DM inventory analysis. The observed DM can be separated into four 
primary components (all in the observer's frame): 


DMs = DMyy + DMyy halo + DMicm + DMhost (3) 


where DMwy is the contribution from the Milky Way’s interstellar 
medium, DM,,), is the contribution from the Milky Way halo, DM, ost 
the contribution from the host galaxy including its halo and any gas 
local to the FRB source, and DM <x, is the contribution from the IGM. 

We use the NE2001 model" to evaluate DMyw = 60 pc cm? (compared 
with 50 pc cm? from the YMW16 model*®) as the mean of a uniform 
distribution with a generous +40% width to conservatively represent 
the uncertainty in DMyy for the direction of FRB 20190520B. For the 
MW halo contribution, we use a uniform distribution from 25 pc cm? 
to 80 pc cm”. Together, the Milky Way disk and halo components yield 
atotal range from 61 pc cm?to164 pc cm? with a mean of 113 pc cm? 
andr.m.s. uncertainty of 17 pc cm”. 

For the IGM, we use the A cold dark matter model (ACDM, where A is 
the cosmological constant) to calculate the mean DM contribution” 
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where the free electron number per baryon in the Universeis y(z) = 7/8 
(which is assumed to be constant for FRB redshifts), the normalized 
matter density is Q,,- 0.315 * 0.007, the dark-energy fraction 
is Q4 c 1- O4, the baryonic fraction is Q, = h ? x (0.02237 x 0.00015), 
the fraction of baryons in the IGM is fi;,, and the Hubble constant is 
Hyzhx100km $12 67.36*0.54 kms ! Mpc! (ref. *). Using values 
for the speed of light c, the gravitational constant G and the proton 
mass M, the resulting expression 
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yields DMig,(0.241) = 248f.., pccm ^. 

Using DM,gy(z) as the mean, we calculate the range of values for a 
given redshift using a log-normal distribution with variance 
CA, = DM. (2DM,, where DM, = 50pccm ?is chosen to provide 
cosmic variance consistent with published simulations. This gives 
Opw ,KO-241) = 111 em pccm ? . The parameters for the log-normal 
distribution are then oj... - { INE + (oy, (2/DMig 2)? BY? and 
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distribution is asymmetric. Equation (5) and the log-normal distribu- 
tion are also used to estimate the median DM,,y, and its inner 68% uncer- 
tainty range shown in Fig. 3. 

To constrain the host-galaxy DM, we combine the MW and IGM 
estimates and their uncertainties with the measured DM averaged 
over all bursts, which isDM,,, = 1,204.7 +4 .Opccm ? The cumulative 
distribution function (CDF) of the log-normal distribution for the 
IGM contribution yields the range DMicy = 195*I P pccm ?. We then 
calculate the posterior distribution of DM, after marginalizing 
over the IGM and Milky Way distributions and using a flat, unin- 
formative prior for DM,,,.. The median and probable ranges are 
calculated from the corresponding CDF. Using fi, = 0.85 for the 
baryon fraction in the IGM'??, we obtain a median value for DM, 
and 68% probable interval DM, = 90372 pccm ? When we vary ficm 
from 0.6 tol, the total range of values for DM, is about 80 pc cm” 
to 350 pc cm ‘and for DM, is about 1,020 pc cm! to 745 pc cm”, 
where we have taken the median values resulting from each value 
of fia, and extended the ranges using (half of the) corresponding 
68% probable ranges for the smallest and largest values of fia. 

For comparison, we also give redshift estimates for different values 
of fi; if a fixed value of DM,,,¢, = 50pccm ^is used (as is often found in 
the literature) along with the above quoted mean value for the MW 
contribution. For fa, ranging from 0.6 to 1, we obtain a range of redshift 
values from 2.2 to 0.9 (including the 68% redshift interval derived from 
the redshift CDF for each value of ficm), much larger than for the redshift 
of 0.241 for the identified host galaxy. 


Constraints on DM,,,. from Ha measurements. The DM of the host 
galaxy is independently estimated from its Ha emission by converting 
the extinction-corrected Ha flux, Fya = (42.0 € 0.5) x 10 “ergs !cm ?, 
to a Ha surface density of 224 + 3 Rayleighs in the source frame at 
z- 0.241, assuming host-galaxy dimensions of 0.5 x 0.5 arcsec (ref. °’). 
These host-galaxy dimensions are only a rough estimate based on the 
size of the Ha-emitting region in the Keck image, but because the image 
is seeing limited, the assumed dimensions may be even smaller, which 
would only serve to increase the Ha surface density in the following 
calculations. The Ha flux and estimated surface density are similar to 
those found for the host galaxy of FRB 20121102A?.. The Ha surface 
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density inthe source frame S(Ha),is related to the Ha surface density 
in the observer frame as S(Ha), = (1+ z)*S(HQ),, and is related to the 
emission measure (EM) in the source frame by 


EM, = 2.75 pc cm °T$?S(Ha), 
S(Ha), | (6) 


z 616-7 pc cm $ x T9? E 


where we have used the redshift z = 0.241 and 7, is the temperature in 
units of 10* K. The EM is related to the DM in the ionized cloudlet 
model" by EM = [C(1** €7)/f]DM7/L, where ( represents cloud-cloud 
variations in the mean density, £ represents the variance of density 
fluctuations in a cloud, fis the filling factor and L is the path length 
through the gas sampled by the FRB. Using this relation, we obtain the 
corresponding source frame DM 


L M floa 
Ç 


1/2 
2 " -3. T045 
DM, = 39243 pc cm °x T; ; kpe d 


xdg (7) 
224t3R| ' 


where we have adopted fiducial values of f= 0.1and = €?- 1, whichare 
typical for the warm ionized medium in the Milky Way. For the path 
length through the gas, we adopt a fiducial value of L =5 kpc, whichis 
based onthe apparentsize ofthe Ha-emitting region in the host galaxy. 
However, given that the optical image is seeing limited, and that the 
orientation ofthe galaxy relative to the FRB line of sight is not known, 
this path length could be as large as 10 kpc or as small as 0.1 kpc. For a 
range of L fromO.1 kpcto10 kpc, we find that DM, could be as small as 
55 pc cm” or as large as 560 pc cm”, for the same fiducial values of T}, 
f, ande’. The estimated DMis also highly sensitive to the temperature: 
for a range of T, from 0.5 to 5, DM, could range from 290 pc cm” to 
810 pc cm”, keeping all other parameters fixed at their fiducial values. 
In the observer's frame, the measured DM contribution of the host 
galaxy is smaller by a factor of 1/(1+ z), yielding a nominal value of 
DMhost,coeft 3165 2pccm ?. for the coefficient in equation (7). The 
quoted errors account for only measurement errors in the Ha flux. 
To match the inferred value of DM,,,, 900 pc cm” requires that the 
three factors in equation (7) involving 7,, L and f combine to a factor 
of about 3, which may easily be explained by a higher temperature or 
the broad range of possible values for L. Regardless, the large Ha EM 
affirms that the FRB DM receives a significant contribution from 
ionized gas in the host galaxy, but itis unclear whether the diffuse, 
Ha-emitting gas can account for the entire host-galaxy contribution 
or whether the FRB's local environment contains significant ionized 
gas content that is not seen in Ha. 


Extragalactic scattering. Scatter broadening manifests as a 
frequency-dependent temporal asymmetry in the burst dynamic spec- 
trum, which is typically modelled as a one-sided exponential pulse 
broadening function convolved with a Gaussian pulse. Some bursts 
from FRB 20190520B have a burst structure that is suggestive of scat- 
tering: leading edges that are aligned across the radio frequency band 
coupled with temporal asymmetries that broaden at lower observing 
frequencies. However, a number of bursts are symmetric across the 
radio frequency band and vary significantly in burst width. When the 
burst S/Nis low, itis also difficultto distinguish scattering from other 
burst substructure, such as the time-frequency drift of intensity is- 
lands? (often referred to as the ‘sad trombone’ effect). Examples of 
bursts with and without evidence of pulse broadening are shown in 
Extended DataFig. 7. 

To measure the mean burst width and scattering timescale of 
FRB 20190520B, we first integrate the dynamic spectrum of each burst 


along the frequency axis and normalize the resulting burst profile. 
Each burst profile is fit with a model composed of a Gaussian com- 
ponent convolved with a one-sided exponential. Out of the 79 fitted 
scattering times, all scattering times with fractional uncertainties 
>50% are excised, leaving a subset of 26 bursts. The scattering times 
of these 26 bursts are then compared with their dynamic spectra 
and one-dimensional burst profiles to verify that the burst temporal 
asymmetry is broadly consistent with a v * frequency scaling. Sup- 
plementary Table 1 shows the scattering times re-scaled to 1.25 GHz 
assuming a v * frequency dependence. The remaining 53 bursts, whose 
scattering times have fractional uncertainties >50%, are then re-fit 
with a Gaussian-only model for the one-dimensional burst profile. 
The full-width at half-maximum burst widths are also shown in Sup- 
plementary Table 1. This provisional approach yields a mean scatter- 
ing time and burst width of 10 + 2.0 ms and 13.5 + 1.2 ms at 1.25 GHz, 
respectively. A more detailed study investigating the burst widths in 
separate frequency subbands is in progress, and corroborates the scat- 
tering interpretation presented here. 

The observed mean scattering time of 10 + 2 msat 1.25 GHz is prob- 
ably too small for the FRB's host galaxy to lie behind our proposed 
host-galaxy association. In this alternative scenario, the host galaxy 
would lie at a redshiftz, > 0.241that depends onthe foreground galaxy's 
contribution to the total DM budget. The FRB would pass through the 
putative intervening galaxy ata redshift z, > 0.241 at animpact param- 
eter of about 4 kpc (based on the observed offset of the FRB localization 
in the optical images). The scattering contribution of the intervening 
galaxy lens is related to its DM contribution by??? 


FG, DME 


T(DM,v,z) = 48.03 usx 1 
( ) p (Le zv 


(8) 


where DM, is the contribution of the intervening galaxy in pc cm? in 
its rest frame, vis the observing frequency in GHz, zis the intervening 
galaxy redshift, and F = (e?/f(/21)!? quantifies the electron density 
fluctuations in the scattering layer for F in units of pc?? km ^, where 
C, e^ and fdescribe the density fluctuation statistics and filling factor, 
lis the outer scale of turbulence and /,is the inner scale. The fluctuation 
parameter È in the Milky Way varies from -10? pc ?? km in the thick 
disk to 210? pc ?? km? near the inner Galaxy, and is typically about 
0.1-1 pc ?? km“? in the thin disk, but it is not generally known for other 
galaxies”. The geometric factor Ge is unity for scattering in a host 
galaxy but it can be very large for an intervening galaxy, in which case 
Gecatt = 2d 4diu/d,, L, where d,, d,, and d,, are the angular diameter dis- 
tances between the source and lens, lens and observer, and source and 
Observer, respectively, and / is the path length through the lens. The 
numerical pre-factor in equation (8) is forallangular diameter distances 
in Gpc and/ in Mpc. 

For the nominal DM contribution implied by the Ha emission, 
DM, = 392 pc cm”, anda path length through the intervening galaxy 
L«5kpc, the implied scattering time is T = 20 s at 1.25 GHz, orders of 
magnitude larger than the observed mean scattering time. For a smaller 
gas temperature T, ~ 0.5, yielding a DM contribution DM, = 290 pc cm”, 
the implied scattering time is still 77 18s at 1.25 GHz, assuming 
F 0.1« pc?? km? and L = 5 kpc. However, there is significant latitude 
inthisestimate, depending on the assumed gas temperature, path length 
sampled by the FRB and the value of È. Although we have assumed that 
the FRB traces the entire DM, implied by the Ha emission, the FRB line 
of sight probably only traces a fraction of the Ha-emitting gas. If DM, is 
assmallas50 pc cm?,thenrcouldbeassmallasO.6 sfor£ 0.1= pc? km ?? 
and/ = 5 kpc. Although the implied scattering time could be two orders 
of magnitude smaller or larger depending onthe combination of F, the 
pathlength/ andthe parameters used to estimate the Ha DM, our fidu- 
cial estimates suggest that the observed scattering time is significantly 
smaller than the scattering expected from an intervening galaxy. 


Rotation measure search from FAST observations. A search for RM 
was performed with the FAST data. The polarization was calibrated 
by correcting for differential gains and phases between the recep- 
tors through separate measurements of a noise diode injected at an 
angle of 45? from the linear receptors. We searched for the RM from 
-30 x 10? rad m?to3.0 x 10? rad m?. Nosignificant peak was found in 
theFaraday spectrum. The observed lack of polarization could be due 
totheintrinsiclow linear polarization, a depolarization process within 
the source or from intra-channel Faraday rotation. 
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and those for the PRS are shown with error bars. These uncertainties include 10 
statistical errors and estimates for systematic errors added in quadrature. 


Extended Data Fig. 1| Positions of the bursts and persistent radio source 
identified with VLA observations. The observation were performed at 1.5,3, 
and 5.5 GHz. They were shown as offsets from the best-fit position of the 
ensemble of bursts, at RA = 16h02m04.272s, dec. = -11? 17’ 17.32" (J2000). The 
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Extended Data Fig. 2| Spectrum of the PRS associated with FRB slightly in order to show the flux density error bars. The inset plot shows the 
20190520B. The bandwidthis split into two subbands for all the observations average flux at each subband. A power-law model fit to these measurements is 


and the corresponding radio flux densities are shown. The observationdatesin shown by the black dashed line and yields a spectral index of-0.41* 0.04 for the 
MJD are shown in the legend. The frequencies for each subband are shifted PRS. 
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Extended Data Fig. 3 | Optical spectrum of the FRB 20190520B host galaxy with FWHM approximately 10 A. The flux scale of the spectrum is not corrected 
obtained at Palomar. The redshift of z= 0.241 was determined usingthe[OIll]- for slit loss or extinction. 
5007À line and Ha line (>50 detections). These two emission lines are narrow, 
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Extended DataFig. 4 | Posterior probability distributions. The probability detected by FAST. b, Results froma fit to all burst waiting times longer than1s. 
distribution are for the shape parameter k and event rater of the Weibull These bursts are at fluences higher than 9 mJy ms. 
distribution. a, Results from a fit to all burst waiting times of FRB 20190520B 
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Extended Data Fig. 5| The fluence-width distribution at1.25 GHzfor FRB 20190520B. The widths and fluences of the 79 bursts detected by FAST are also 
showninSupplementary Table 1. 
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Extended Data Fig. 6| The distribution of burst energies for FRB respectively. The 90% completeness threshold uses 0.023 Jy ms, which is the 
20190520B. The bursts are from FAST observations. The red line and black simulation result from ref. ””. 
dashed line are the log-normal fit and 90% completeness threshold, 
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Extended Data Fig. 7 | Dynamic spectra of bursts with and without signal-to-noise ratio (S/N). For P31, there is a potential combination of 
scattering. a, Frequency-time dynamic spectra of bursts with significant scattering, time-frequency drift, and multiple unresolved burst components; 
evidence of pulse broadening, along with1D burst profiles that have been inthis case we report a scattering time that should be interpreted as an upper 
averaged across the entire frequency band. b, Examples of bursts without limit. Bursts with scattering time constraints are shown in Supplementary 
significant evidence of pulse broadening. White lines indicate excised radio Table 1. 


frequency interference. All ID burst profiles are shownin units of the 
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Frequency refers to the centre of the observing band and duration represents the total length of on-source observations. Note that the realfast system was no 


Extended Data Table 1 | VLA observations of the PRS and number of bursts detected in each observation 
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Extended Data Table 2 | Localized positions of the PRS from 1.5, 3, 5.5GHz VLA deep images 


Frequency R.A. Decl. Statistical Error Systematic Offsets 
(GHz) R.A. ("^") Decl. (^) R.A. (") Decl. () 


L5 16h02m04.2543s  —11dl7ml7.4146s 0.0386 0.0713 0.0373+0.2460 0.0144+0.3035 
3 16h02m04.2611s  —11d17ml17.3869s 0.0176 0.0168 —0.0696+0.2772 0.0382+0.1552 
3.5 16h02m04.2618s  —11dl7ml7.3375s 0.0075 0.00905 -—0.0947+0.1161 0.1069+0.0554 


The first column shows the observing central frequency. The second and third columns report the coordinates of the PRS from the deep images. The remaining columns show the 10 statistical 
errors and the cross-matched systematic offsets in right ascension and declination. 


Extended Data Table 3 | Localized positions of the VLA bursts 


Band S/N R.A. Decl. Statistical Error Systematic Error 
R.A.(") Decl. (^) R.A.(") Decl. (”) 
L1 16 16h02m04.2742s  -11d17m17.0535s 0.124 0.208 0.2460 0.3035 
L2 27  16h02m04.2847s  -11dl7ml7.1912s 0.060 0.145 0.2460 0.3035 
L3 15  16h02m04.2691s -11d17m17.1299s 0.070 0.139 0.2460 0.3035 
S1 19 16h02m04.2567s -—11d17m17.4093s 0.044 0.080 0.2772 0.1552 
S2 17  16h02m04.2449s .-11dl7ml7.1874s 0.153 0.173 0.2772 0.1552 
S3 15  16h02m04.2602s .—11dl7ml17.5078s 0.103 0.176 0.2772 0.1552 
S4 10  16h02m04.2661s .-—11dl7ml7.3666s . 0.333 0.247 0.2772 0.1552 
S5 7 16h02m04.2500s -11dl7ml17.9337s 0.214 0.269 0.2772 0.1552 
Cl 19  16h02m04.2734s -11dl7ml7.3078s 0.041 0.071 0.1161 0.0554 


The label in the first column shows the frequency band of observation, followed by the burst number. DM represents the S/N maximizing DM obtained using offline refinement of the bursts and 
S/N reports the maximum image S/N. The remaining columns show the burst positions (and errors) obtained using CASA calibration and image fitting. The systematic errors reported here were 
estimated using deep radio images generated using all the observations at the respective frequency band. 
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Extended Data Table 4 | Properties of the VLA bursts 


Name MJD DM Flux density (Jy) Width (ms) Frequency (GHz) 
L1 59051.0476715 1209.60 0.113+0.006 20 1.33-1.52 
L2 . 59053.1045803 1227.15 0.275+0.010 20 1.39-1.52 
L3 59053.0973272 1236.60 0.189+0.009 10 1.78-1.97 
S1 59091.9604318 1291.30 0.075+0.004 20 3.12-3.50 
S2 59091.9560792 1193.40 0.100+0.009 20 2.63-2.88 
S3 59091.9635008 1222.75 0.093+0.007 20 2.50-2.76 
S4 59091.9435287 1216.80 0.037+0.005 10 2.88-3.24 
S5 59091.9668476 1276.50 0.022+0.003 10 2.50-3.50 
Cl 59107.9272138 1267.50 0.068+0.004 10 4.49—4.88 


The label in the first column shows the frequency band of observation, followed by the burst number. MJD is referenced to the solar system barycentre (in the TDB scale) and corrected for 
dispersion delay at infinite frequency. Width values here should be considered as upper limits. Frequency represents the range of frequencies in which the burst signal is prominently seen, and 
was visually determined. The DM that maximizes the S/N of the respective burst is given in the third column. The apparent DMs may reflect time-frequency structure, as opposed to bona fide 
DM variations. 


